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^ ■ We calculate bulk thermodynamic properties, such as the pressure, energy density, and entropy, 

^ , in SU(4) and SU(8) lattice gauge theories, for the range of temperatures T < 2.{)Tc and T < l.GTc 

I respectively. We find that the = 4, 8 results are very close to each other, and to what one finds 

I in SU(3), and are far from the asymptotic free-gas value. We conclude that any explanation of the 

■ high-T pressure (or entropy) deficit must be such as to survive the N ^ oo limit. We give some 
examples of this constraint in action and comment on what this implies for the relevance of gravity 

, duals. 
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I. INTRODUCTION 



The thermodynamic properties of Quantum Chromodynamics (QCD), besides being of 
fundamental interest, are currently at the centre of intense experimental research. One of 
the most interesting phenomena has to do with the range of temperatures, T, above the 
phase transition (or crossover) at T = Tc, where the theory deconfines and chiral symmetry 
is restored. Traditionally, the description of this transition assumed that the hadronic phase 
gives way to a plasma, whose physical degrees of freedom are weakly interacting quarks 
and gluons. Recent experimental results have, however, challenged this 'simple' picture (for 
example see jlj and references therein), and point to a picture of the 'plasma' as a very good 
fluid in the accessible range of T above T^. In fact, numerical lattice results had already 
demonstrated the inadequacy of the simple quark-gluon plasma picture some time ago. Such 
lattice calculations, both for the pure gauge case [2] and with different kinds of fermions j3| , 
found a large deficit in the pressure and entropy as compared to the Stephan-Boltzmann 
predictions for a free gluon gas (for pure glue), which remained at the level of more than 
10% even at temperatures as high as T ~ 4Tc. Further evidence that points in the same 
direction is the survival of hadronic states above T^, as seen in recent lattice simulations (for 
example see 0] and references therein). 

These lattice calculations, and more recent experimental observations, have attracted 
considerable attention (see e.g. i5] for a review). Approaches have ranged from modeling 
the system in terms of noninteracting quasi-particles with the quantum numbers of quarks 
and gluons but with temperature dependent masses 0,01, to using higher order perturbation 
theory (restricted by infrared divergences), sometimes including nonperturbative contribu- 
tions on the dimensionally reduced 3D Euclidean lattice 0|, large re-summations (e.g. ^ 
and references therein), or, more recently, a description |lO| in terms of a large number of 
loosely bound states that survive deconfinement and come in various representations of the 
gauge and flavor groups, and where one can use for example the lattice masses measured in 
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In this paper we ask whether this pressure (and entropy) deficit is a dynamical feature 
not just of SU(3) but of all SU(A^) gauge theories - and in particular whether it survives 
the N —y oo limit. In this limit the theory becomes considerably simpler, although not 
(yet) analytically soluble, and so what happens there should strongly constrain the possible 
dynamics underlying the phenomenon. For example, in that limit supersymmetric SU(A^) 
gauge theories become dual to weakly coupled gravity models, and in that context we recall 
the frequently mentioned prediction 2|, that the pressure in the strong-coupling limit of 



the J\f = 4 and N = oo supersymmetric gauge theory is 3/4 of its Stephan-Boltzmann value, 
which is similar to the deficit, referred to above, that one finds in the non-supersymmetric 
case. 

To address this question we calculate the pressure for T < 2Tc in S'f/(4), and SU{8) 
lattice gauge theories and compare the results to similar SU(3) calculations available in the 
literature (which we supplement where it is useful to do so). Recent calculations of various 
properties of SU(A^) gauge theories [l3| have demonstrated that SU(8) is in fact very close 
to SU(cxo) for most purposes and have provided information on the location, /3c, of the 
deconfining transition for various Lf and N [McilSj]. Thus our calculations should provide 
us with an accurate picture of what happens to the pressure at = oo. 

In the next Section we summarise the lattice setup, the relevant thermodynamics, and 
provide numerical checks that our system is large and homogeneous enough for our thermo- 
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dynamic relations to be appropriate. We then present our results for the pressure, entropy 
and related quantities. We discuss the implications of our findings in the concluding section. 



II. LATTICE SET UP AND METHODOLOGY 

The theory is defined on a discretised periodic Euclidean four dimensional space-time with 
X Lt sites. Here Ls,t is the lattice extent in the spatial and Euclidean time directions. 
The partition function 

Z{T, V)=J: exp [-^] = exp {-^} = exp |-^| (2.1) 

defines the free energy F and the free energy density, /, and can be expressed as a Euclidean 
path integral 

Z{T,V) = J Df/exp(-/35w). (2.2) 

Here T = {aLt)~^ is the temperature and V = (aL^)^ is the spatial volume. When we 
change /?, so as to change the lattice spacing a{f3), we change both T and V, if Lg and Lt 
are kept fixed. In the large- limit, the 't Hooft coupling A = g'^N is kept fixed, and so we 
must scale (3 = 2N'^/X oc A^^ in order to keep the lattice spacing fixed in that limit. We use 
the standard Wilson action S'w given by 



p 



1 - ^ReTr Up 
N 



(2.3) 



Here P is a lattice plaquette index, and Up is the plaquette variable obtained by multiplying 
link variables along the circumference of a fundamental plaquette. We perform Monte-Carlo 
simulations, using the Kennedy- Pendelton heat bath algorithm for the link updates, followed 
by five over-relaxations of all the SU{2) subgroups of SU{N). 



A. The method used 

In lattice calculations of bulk thermodynamics, one can choose to use either the "integral" 
method (e.g. j3]) or the "differential" method (e.g. or a new variant jl3]) or one can 
attempt a direct evaluation of the density of states (e.g. [3|). We choose to use the first 
of these methods since the numerical price involved in using larger values of N drives us 
to smaller Lt, which means that the lattice spacing is too coarse (about O.lSfm) for the 
differential method. We have performed preliminary checks for the applicability of the 
Wang-Landau algorithm [l^ for the evaluation of the density of states in the SU{8) gauge 
theory, but found it numerically too costly for the present work. 

The properties we will concentrate on are the pressure p, the energy density per unit 
volume e, and the entropy S, as a function of temperature. These are given by 

P = T-^ log Z{T, V) = ^ log Z{T, V) = -/, (2.4) 

6 = ^AiogZ(T,\/), (2.5) 

^ = = i±^. (2.6) 

V T T ^ ^ 
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where the second equahty in the first and last hnes follows if the system is large and homo- 
geneous, i.e. if V is large enough. In addition it is useful to consider the quantity 



.5 d P 



A^e-3p = T'—jr, (2.7) 

which vanishes for an ideal gluon plasma. Again the second equality requires a large enough 
V. To calculate the pressure at temperature T in a volume V with lattice cut-off a(/5), we 
express logZ in the integral form 

p{T) = -logZ{T,V)= / d(3'—l^ (2.8) 

(There is in general an integration constant, but it will disappear when we regularise the 
pressure later on in this section.) This integral form is useful because it is easy to see from 
Eqs. fl2.2l2.;^|) that 

= -{Sw) = N,{u,) (2.9) 

where Np = GLfL^ is the total number of plaquettes and Up = ReTr Up/N. So the pressure 
can be obtained by simply integrating the average plaquette over j3. This pressure has been 
defined relative to that of the unphysical 'empty' vacuum and will therefore be ultraviolet 
divergent in the continuum limit. To remove this divergence we need to define the pressure 
relative to that of a more physical system. We shall follow convention and subtract from 
p(T) its value at T = 0, calculated with the same value of the cut-off a(/5). Thus our pressure 
will be defined with respect to its T = value. Doing so we obtain from Eq. ()2.91 12. 8|) 



a'\p{T) -p{0)] = 6 f d[3\{up)T - {up)o). (2.10) 



where (Mp)o is calculated on some lattice which is large enough for it to be effectively at 
T = 0. We replace p{T) — p(0) p{T), where from now on it is understood that p{T) is 
defined relative to its value at T = 0, and we use T = {aLt)~^ to rewrite Eq. ()2.10j) as 

^^^^ - 6Lt fd(3\{up)T - {up)o). (2.11) 



We remark that when our L^Lt lattice is in the confining phase, then (up) is essentially 
independent of Lt and takes the same value as on a lattice (see below). This should 
become exact as ^ cxd but is accurate enough even for SU(3). Thus as long as we choose 
Po in Eq. ()2.1ip such that a{Po)Lt > l/T^ then the integration constant, referred to earlier, 
will cancel. 

Finally, we evaluate A in Eq. ()2.7|) as follows: 

A/T^ = (2.12) 
^ P (2.13) 



d\ogTdl3T^ 



QLtiiupmo - (upmr) X ^i^l^- (2.14) 
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To evaluate d\og(a(p))/dp we can use calculations of the string tension, a, in lattice units. 
For example, in [20| the calculated values of ay/ a are interpolated in P for various N and 
one can take the derivative of the interpolated form to use in Eq. ()2.14|) . One could equally 
well use the calculated mass gap or the deconfining temperature. All these choices will give 
the same result up to modest O(a^) differences. 



B. Average plaquette 

We see from the above that what we need to do is to calculate average plaquettes closely 
enough in [3 so as to be able to perform the numerical integration in j3. And we need the 
average plaquettes not only on the LtL^ lattice but also on a reference 'T = 0' lattice at 
each value of /3. However we mostly need values for /3 > Pc, where a{(3c)Lt = l/Tc, since 
p{T) - p{0) ^ once T < T^. 

We performed calculations in SU{4) on 16^5 lattices and in SU{8) on 8^5 lattices for a 
range of (3 values corresponding to T/T^ e [0.89, 1.98] for SU{A), and to T/T^ e [0.97, 1.57] 
for SU{8). Since we use Lt = 5, while the data for SU{3) in [2| is for Lf = 4,6,8, we 
also performed simulations for SU{3) on 20'^5 lattices with T/T^ G [1,2]. The results are 
presented in Tables HV lIIII 

In addition to the finite T calculations we have performed 'T = 0' calculations on 20^ 
lattices for SU(3), and on 16^ lattices for SU(4). These have the advantage of being on the 
same spatial volumes as the corresponding finite T calculations, and we know from previous 
calculations [2l[|23] that, for the range of a(/?) involved, these volumes are large enough to 
be, effectively, at zero T. For SU(8) however, using 8^ lattices would not be adequate for the 
largest /^-values, as we will see below. (The same is not true for the finite T calculation on 
8'^5 lattices where it is 1/aT that sets the scale for finite volume corrections.) We therefore 
take instead the SU(8) calculations on larger lattices in [i^l, and interpolate between the 
values of P used there, to obtain average plaquettes at the values of /5 we require. To perform 
this interpolation we fit with the ansatz 

(upUP) = {u,)nP) + ^]^(«v^)' + + c,g'', (2.15) 



where {up)q'^'{P) is the lattice perturbative result to 0{g^) from j23| and = 8. Our best 
fit has x^/dof = 0.93 with dof = 2, and the best fit parameters are C4 = —6.92, C5 = 26.15, 
and a gluon condensate of = 0.72. 

For the scaling of the lattice spacing with /5, needed in Eq. fl2.15|l and Eq. ()2.14|1 and in 
the temperature scale, we used the interpolation of ay/a as a function of /3, as given in 
^. For the temperature scale we need in addition to locate the value of P that corresponds 
to T = Tc for the relevant value of Lt, and for this we have used the values in 3, 2^- In the 



case of SU{3) we compared the resulting T/Tc{P) with that of [2] where the physical scale 
was set by T^. We find that the two functions lie on top of each other for = 6. This is 
consistent with the fact that the SU (3) value of Tc/ yfa for Lj = 5, 6 are the same within one 



^ This is excluding the first three /3 values in the case of S'C/(4), which are outside the interpolation regime 
of [i^. In that case we have performed a new interpolation fit to include these points as well. This 
gave the string tensions ayfa = 0.3739(15), 0.3440(10), 0.3336(10) and the derivatives —dloga/d(3 = 
1.83(7), 1.55(7), 1.48(5) for /3 = 10.55, 10.60, 10.62. 
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sigma [13 • This is true for SU{8) as well, where the value of T^/ y/a for a = l/(5Tc), and 
a = l/(8Tc), are the same within one sigma j23], and we find no point to perform similar 
comparisons there. For SU^i) the value ofTl/y/a at a = l/(5Tc), l/(6Tc) is ~ 5, and ~ 3.7 
sigma away from the value at a = l/(8Tc) |j20|, which may suggest that in this case T/Tc{(3) 
at values of (3 that correspond to T ~ 8/5Tc will be smaller when fixing the physical scale 
with Tc rather than with the string tension. Nevertheless the shift between the two is at the 
level of ~ 2%, and will not change the results presented here. In addition, to fix T/Tc(/3) by 
fixing Tc, requires a larger scale calculation of PdLt, Lg) that will include evaluation of finite 
volume corrections, similar to what was done for = 5 in 3- In view of the small shifts 



and the high calculational price, we shall ignore this potential ambiguity in this paper. 



TABLE I: Statistics and results of the Monte-Carlo simulations for 5'C/(4). 





T > 


T = 




ST 


(lattice sweeps) x 10 ^ 


So 


(lattice sweeps) x 10 


10.55 


0.537478(84) 


10 


0.537487(81) 


5 


iU.DU 


0.543862(58) 


20 


0.543797(25) 


15 


10.62 


0.546212(64) 


10 


0.546068(33) 


10 


10.64 


0.550279(70) 


10 


0.548208(16) 


20 


10.68 


0.554213(32) 


20 


0.552177(16) 


20 


10.72 


0.557649(30) 


20 


0.555861(14) 


20 


10.75 


0.560051(27) 


20 


0.558462(13) 


20 


10.80 


0.563923(32) 


20 


0.562587(16) 


20 


10.85 


0.567592(24) 


20 


0.566453(17) 


20 


10.90 


0.571107(17) 


20 


0.570118(16) 


20 


11.00 


0.577707(17) 


20 


0.576981(11) 


20 


11.02 


0.578985(18) 


20 


0.578279(11) 


20 


11.10 


0.583911(20) 


20 


0.583352(12) 


20 


11.30 


0.595398(13) 


20 


0.595039(10) 


20 



C. Finite volume effects 

For = 4, 8, one is able to use lattice volumes much smaller than what one needs for 
SU {2>) \M. That this is so for the deconfinement transition, has been explicitly demonstrated 
in [l^ |20|, and is theoretically expected, much more generally, as — > 00. The main 
remaining concern has to do with tunnelling between confined and deconfined phases near T^. 
When ^ 00 tunnelling occurs only aX (3 = (3c (in a calculation of sufficient statistics) and 
the system is in the appropriate pure phase for T < and for T > T^. On a finite volume, 
where this is no longer true, one minimises finite- \^ corrections by calculating the average 
plaquettes only in field configurations that are confining, for T < Tc, or deconfining, for 
T > Tc. This ensures that the system is as close as possible to being 'large and homogeneous' 
as is required in the derivations of this Section. Because the latent heat grows oc iV^ 
the region ST around Tc in which there is significant tunnelling shrinks as ST oc for 
a given V. Hence we can reduce as iV increases without increasing the ambiguity of 
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TABLE II: Statistics and results of the Monte-Carlo simulations for SU{8). 



T > 


T = 


P 


ST 


Ls 


(lattice sweeps) x 10 


/3 


So 


(lattice sweeps) x 10 


43.90 


0.525330(80) 


14 


5 


43.85 


0.523819(37) 


> 20 


43.93 


0.526873(79) 


8 


19.5 


44 


0.528788(18) 


> 20 


44.00 


0.531307(50) 


10 


> 20 


44.35 


0.538491(13) 


> 20 


44.10 


0.534164(34) 


12 


7 


44.85 


0.549794(9) 


> 20 


44.20 


0.536650(70) 


14 


5 


45.7 


0.565708(4) 


> 20 


44.30 


0.539181(30) 


8 


20 








44.45 


0.542629(38) 


8 


30 








44.60 


0.545812(35) 


8 


20 








44.80 


0.549968(37) 


8 


30 








45.00 


0.553926(38) 


8 


20 








45.50 


0.562992(28) 


12 


10 









TABLE III: Statistics and results of the Monte-Carlo simulations for SU{3). 



(3 


T > 


T = 




ST 


(lattice sweeps) X 10 ^ 


So 


(lattice sweeps) x 10 


5.800 


0.568664(100) 


10 


0.567667(29) 


11 


5.805 


0.569688(153) 


20 


0.568438(23) 


11 


5.810 


0.570624(55) 


10 


0.569218(18) 


11 


5.815 


0.571297(81) 


10 


0.569996(26) 


11 


5.820 


0.572205(78) 


10 


0.570788(16) 


11 


5.900 


0.583058(38) 


10 


0.581854(20) 


11 


6.150 


0.609377(27) 


10 


0.608971(8) 


11 


6.200 


0.613966(31) 


10 


0.613628(13) 


11 



the calculation. For SU(3), where the phase transition is only weakly first order, frequent 
tunneling occurs in the vicinity of Tc in the volume we use, and it is not practical to attempt 
to separate phases. This will smear the apparent variation of the pressure across in the 
case of SU(3). 

We now turn to a more detailed discussion of finite volume effects. If ^ is the longest 
correlation length, in lattice units, in a volume of length L, then finite volume effects will 
be negligible if <^ L. In addition finite volume corrections will be suppressed as ^ oo. 
In our particular context, ^ is given by the inverse mass of the lightest (non- vacuum) state 
that couples to the loop that winds around the temporal torus. In both the confined and 
deconfined phases, these masses decrease as T ^ Tc. Therefore the largest length scale is 
set by the masses at T = Tc. As increases these masses increase towards their limits, 
with l/N"^ corrections that are quite large (2ol |. 
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1. The deconfined phase 



In the deconfined phase, on an x 5 lattice at T = T^, the value of ^ is about 12.5 
lattice spacings for 577(3), while it is about 5.2, and 2.4 lattice spacings for SU{4:), and 
SU{8) respectively This suggests that our choice of Ls = 16 for SU(4) and Ls = 8 for 
SU(8) should be adequate. In addition it is known from calculations of Tc j^EHE^ ^^at on 
such lattices the tunnelling is sufficiently rare that even at T = Tc one can reliably categorise 
field configurations as confined or deconfined and hence calculate the average plaquette in 
just the deconfined phase if one so wishes. For our supplementary SU(3) calculations we use 
Ls = 20 which is much smaller in units of ^. In practice this means that in this case we are 
unable to separate phases at T ~ T^. 

To explicitly confirm our control of finite volume effects, we have compared the SU(8) 
value of {up{(3)) as measured in the deconfined phase of the our 8^ x 5 lattice with other 

X 5 results from other studies [l^. As summarised in Table HV| the results are consistent 
at the 2cT level. 



TABLE IV: Finite volume effects for plaquette average in the deconfined phase on a Lf = 5 lattice, 
for 5f/(8). 





Ls = ^ 


Ls = 10 


Ls = 12 


Ls = 14 


43.95 




0.529788(100) 


0.529944(65) 




44.00 


0.531343(45) 


0.531307(50) 






44.10 


0.534219(54) 




0.534164(34) 




44.20 


0.536714(33) 




0.536689(54) 


0.536650(70) 


44.25 






0.537954(60) 


0.537850(100) 


44.30 


0.539181(29) 






0.539220(100) 


45.50 


0.563093(41) 




0.562992(28) 





2. The confined phase 

As we remarked above (see below for explicit evidence) we have (mp)t — (wp)o in the 
confined phase and so the contribution in Eq. ()2.1H1 of the range of (3 where the finite T 
system is confined is very small. Nonetheless, we include an integration over that range for 
completeness and so we need to discuss possible finite V corrections for this well. 

In the confined phase, on an x 5 lattice at T = T~, the value of ^ is about 9.5 lattice 
spacings for SU{?>), but drops to about 5 and 3.5 for SU{4) and SU{8) respectively |2Q|]. 
This leaves our choice of Lg still reasonable for S'[/(4) but somewhat worse for SU(8). In 
Table |3 we provide a finite volume check for the latter case that proves reassuring. 

Finally we return to our earlier comment that for the 'T=0' L^ lattice calculations, a size 
L = 8 in SU(8) would not be large enough. This is demonstrated, for our largest /5-value, 
in Table IVlj where we also present the value of Lt x T /Tc{(3) (in our Lt = 5 calculations). 
In the confined lattice, finite volume effects will be suppressed when the latter is much 
smaller than L^. Clearly for j3 = 45.70, and Ls = 8, this is not the case. 

By contrast, for SU{4) the finite volume effects seems not to be large on the 16^ lattice 
as we checked for our largest value of (3 = 11.30. There the value of the plaquette on a 20"^ 
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lattice is 0.595014(4) [21|, which is consistent within ~ 2.3 sigma with the value presented 
in Table m This is in spite of the fact that for this coupling Lt x T/T^ = 10, and is not so 
much smaller than = 16. 



TABLE V: Finite volume effects for plaquette average in the confined phase on a = 5 lattice, 
for SC/(8). 



/3 


= 8 


Ls = 10 


Ls = 12 


Ls = 14 


43.90 


0.525750(87) 




0.525613(54) 


0.525425(90) 


43.95 




0.527240(34) 


0.527275(48) 


0.527280(50) 


44.00 






0.528867(33) 


0.528810(50) 


44.10 






0.531880(45) 


0.531900(60) 



TABLE VI: Finite volume effects for plaquette average in the confined phase on a lattice, for 
SU{8). The last column is for Lt = 5. 



13 


Ls = 8 


Ls = 10 


Ls = 16 


Lt X T/T, 


44.00 


0.528876(39) 


0.528788(18) 




5.05 


45.70 


0.566089(23) 




0.565708(4) 


8.20 



III. RESULTS 

To obtain the pressure from the values of the average plaquette presented in Tables HV lIIII 
we need to perform the integration in Eq. ()2.1H) . which we do by numerical trapezoids. 
We have already remarked that the contribution to the pressure from the confined phase is 
negligible. In Table IVIII we provide some accurate evidence for this. We show the values 
of the average plaquette on L^ lattices, corresponding to T ~ 0, as well as the values on 
Lg5 lattices at T ~ T^, with the latter obtained separately in the confined and deconfined 
phases. (These volumes are large enough for there to be no tunnelling, or even attempted 
tunnelling, within our available statistics.) We see that for both SU{4) and SU{8) there is 
no visible difference between the plaquette at T = and T = Tc in the confined phase at, 
say, the 2a level. Any difference, (and there obviously must be some difference) is clearly 
negligible when compared to the difference between the confined and deconfined phases at 
(and above) T^. 

In presenting our results for the pressure, we shall normalize to the lattice Stephan- 
Boltzmann result given by 



Here Rj includes the effects of discretization errors in the integral method 2^ 2^- For large 
values of Lt, and an infinite volume, it is given by 

ML.)^l^l(^)\l(^)\0(r)\ (3.2) 



21 \LtJ 21 \LtJ \Lt 
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TABLE VII: The plaquette average in the confined phase, C, at T ~ Tc compared to the T = 
value and to the value in the deconfined phase, D. For S'C/(4) and SU{8). 



f3 




lattice 




phase 


T 


10.635 


4 


32^5 


0.549563(33) 


D 


Tc 






32^5 


0.547689(11) 


C 


Tc 






iU 


U.04 ( D4U(^Z ( j 




U 


43.965 


8 


12^5 


0.530352(23) 


D 


Tc 






12^5 


0.527725(27) 


C 


Tc 






10^ 


0.527648(24) 


C 






Since some values of Lt discussed in this work are not very large, we shall use the full 
correction, which includes higher orders in 1/^, instead of Eq. ()3.2j) . This was calculated 
numerically for the infinite volume limit in |25| for Lt = 4, 6, 8, and we supplement this 
calculation, with the same numerical routines [26[, for other values of L^. A summary of 
Ri{Lt) in the infinite volume limit is given in Table IVTTTl 



TABLE VIII: The lattice discretisation errors correction factor Ri{Lt) in the infinite volume limit. 



Lt = 2 


it = 3 


Lt = A 


Lt = b 


Lt = Q 


it = 8 


2.04526(4) 


1.6913(2) 


1.3778(1) 


1.2129(6) 


1.1323(1) 


1.0659(1) 



We find that the full correction for Lt = 5 is a ~ 21% effect, which, without this normal- 
isation, might obscure the physical effects that we are interested in. This is an appropriate 
normalisation because we expect Eq. (jH.lj) to provide the T oo limit of p/T^. The same 
applies to the internal energy density, since e — *• 3p as T ^ cxd, and so when presenting our 
results for e/T^ we normalise it with the expression in Eq. p.ip . For similar reasons we shall 
use the same normalisation when presenting our results for the entropy. For A/T^ it is less 
clear what normalisation one should use since A = e — 3p^0asT^oo, but for ease of 
comparison we shall once again normalise using Eq. ()3.1|1 . 

To facilitate the comparison of our results with earlier work on SU(3) [2J, which was 
done for Lt = 4,6,8, we have performed SU{3) simulations with Lt = 5. The spatial size 
is Ls = 20 which should be sufficiently large in the light of our above discussion of finite 
volume effects (and we note that it satisfies an empirical rule that one needs L^/Lj > 4 (27l|). 

We present our = 4 and A^ = 8 results for p/T^ in Fig. ^ We also show there our 
calculations of the SU{3) pressure for Lt = 5, as well as the Lt = 6 calculations from 0. 
Although our errors on the SU{3) pressure are probably underestimated, since the mesh in 
P is quite coarse, nonetheless one can clearly infer that the pressure in the S'f/(4) and SU (8) 
cases is remarkably close to that in SU{3) and hence that the well-known pressure deficit 
observed in SU(3) is in fact a property of the large- A^ planar theory. 

In Fig. 121 we present our results for A/T^ as calculated from Eq. ()2.14|1 . This quantity 
can be considered as a measure of the interaction and non-conformality of the theory, since 
it is identically zero both for the noninteracting Stephan-Boltzmann case, and for the A/" = 4 
super symmetric SU{N) gauge theory. As remarked above, we normalise with the expression 
in Eq. ()3.1|) . We also note that in this case there are no errors from a numerical integration, 
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FIG. 1: The pressure, normalized to the lattice Stephan-Boltzmann pressure, including the full 
discretization errors. The symbol's vertical sizes are representing the largest error bars (which are 
received for the highest temperature). The solid line is for SU{3) and Lt = 6 from 

and this enables a fair comparison with the SU{3) data of 0. Comparing the results for 
different N we see that, just as for the pressure, the results for all these gauge theories are 
very similar. 

To see what is the behaviour of A/T^ at even higher temperatures, we use the plaquette 
averages on lattices with Lf = 2, 3, 4, 5, that have been calculated at fixed couplings which 
correspond to T ~ Tc for Lf = 5 [2i|. We present the results in Table HXl For the evaluation 
of A one needs d\oga/d[3 which we present in the table as well. 

TABLE IX: Plaquette average in the deconfined phase for lattice with fixed coupling, at different 
values of L^, and with (5 that corresponds to roughly the deconfining temperature at Lt = 5: 
13 = 5.800, 10.635, 44.00 for = 3, 4, 8. The data for = 5 are obtained for L = 64, 32, 10 for 
A = 3, 4, 8 (for A = 3, 5{up) is the difference between the plaquette as calculated within separate 
confined and deconfined sequences of field configurations). 



N 


L3 X 5 


8^ X 4 


8^ X 3 


8^ X 2 


10^ 


—d\oga/d[3 


3 


5{up) = 0.00080(5) 


0.570987(37) 


0.573311(34) 


0.578121(27) 


0.567642(29) 


2.075(17) 


4 


0.549563(33) 


0.551604(33) 


0.554047(27) 


0.559163(24) 


0.547640(27) 


1.440(23) 


8 


0.531202(92) 


0.533066(25) 


0.535991(24) 


0.541518(17) 


0.528788(18) 


0.384(20) 



In such calculations where one varies T by varying Lt, the lattice spacing varies as a = 
X 1/T when expressed in units of the relevant temperature scale, and so lattice spacing 
corrections will vary with T. 

The resulting values of A in the case of SU(3) are plotted in Fig. El where they are 
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FIG. 2: Results for A{T)/T^ = T^^, normalized by the same coefficient as we normalize the 
pressure. The sohd hne is for 5C/(3) and Lt = 6 from 

compared to the results obtained from calculations where one varies T by varying (3 at fixed 
Lf. These calculations include ours for Lt = 5 and those of [3] for = 4, 6. 

As we see from Fig.lHlour Lt = 5 SU (3) results do in fact lie between the Li = 4, 6 results 
of as one would expect. We observe that the T dependence is very similar in all cases, 
and that the remaining Lt dependence appears to be much the same for the different kinds 
of calculation. This gives us confidence that performing calculations where we vary T by 
varying Lt at fixed (3 does not introduce any unanticipated and important systematic errors. 

Having performed this check, we compare in Fig. |3 our results for A in the range < 
T < 2.5Tc that corresponds to 5 > > 2. This comparison confirms what we observed in 
Fig. 121 over a smaller range of T: A is very similar for all the values of (except very close 
to Tc), implying that this is also a property of the N = oo planar limit. 

Finally we present in Fig. our results for the normalized energy density e = A + 3p, and 
the entropy per unit volume s = {e+p)/T. The lines are the SU{3) result of with Lt = 6. 
Again we see very little dependence on the gauge group, implying very similar curves for 
A^ = oo. 



IV. SUMMARY AND DISCUSSION 



In this work we have analyzed numerically the bulk thermodynamics of SU(4) and SU(8) 
gauge theories. We found that the pressure, when normalized to the Stephan-Boltzmann 
lattice pressure, is practically the same as for SU(3), in the range Tc < T < l.GTc that we 
analyze. We found the same to be the case for the internal energy and entropy, as well as for 
the quantity A = e — 3p (where we were able to explore temperatures up to T ~ 2.5Tc). All 
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FIG. 3: Results for MT)/T^ = T^^, normalized to the free-gas result. The lines are for SU{3) 
and Lt = 4,6 from Red triangles correspond to Lt = 5, and changing (3, while blue circles 
correspond to changing Lt and keeping a fixed /? = 5.800. 

this implies that the dynamics that drives the deconfined system far from its noninteracting 
gluon plasma limit, must remain equally important in the N = oo planar theory. This is 
encouraging since that limit is simpler to approach analytically, in particular using gravity 
duals. 

Our results have been (mostly) obtained for lattice spacings a = 1/{5T) and it would be 
useful to perform a larger scale calculation that allows us to perform an explicit continuum 
extrapolation. However past SU(3) calculations of the pressure, and calculations in SU(A^) 
of various physical quantities, strongly suggest that our choice of a already provides us with 
a reliable preview of what such a more complete calculation would produce. 

Our results imply that any explanation of the QCD pressure deficit must survive the 
large- limit, and so should not be driven by special features particular to SU(3). This 
can provide a strong constraint on such explanations. For example, in approaches based 
on higher order perturbation theory, it tells us that the important contributions must be 
planar. In models focussing on resonances and bound states, it must be that the dominant 
states are coloured, since the contribution of colour singlets will vanish as — oo. Models 
using 'quasi-particles' should place these in colour representations that do not exclude their 
presence at = cxd, and in fact give them T-dependent properties which depend weakly on 
N. Also, topological fluctuations should play no role in this deficit since the evidence is that 
there are no topological fluctuations of any size in the deconfined phase at large- A^ [28, 2_9J . 

Finally, we emphasize that our conclusion that the SU(3) pressure and entropy deficits 
are features of the large-A^ gauge theory, means that these 'observable' phenomena can, 
in principle, be addressed using AdS/CFT gravity duals. Indeed it is precisely where the 
deficit is large that the coupling must be strong and this is also precisely where, at large 
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FIG. 4: Results for A{T)/T^ = for N = 3,4,8, by fixing (3 = (3c{Lt = 5), while changing 

Li = 2,3,4,5. 

iV, such dualities can be established. As has been frequently emphasized (see for example 
[id Il7| ) the deficit in the normalized entropy is not far from the value of s/sfree-gas = 3/4 
given by the AdS/CFT prediction. In this paper we have found that large- gauge theories 
show the same behaviour, as we see in Fig. El where, for the entropy, the horizontal line 
-Snormaiized/T"^ = 3 would corrcspoud to s/sfree-gas = 3/4. Our rcsults cau therefore serve as a 
bridge between the AdS/CFT approach to large- A^ and the observable world of QCD. 
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